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Abstract 

In this note we compare the latest 1.04 fb~^ LHC searches for squarks and gluinos from 
jets and missing transverse momentum (MET) with constraints arising from the LEP Higgs 
bound. For General Gauge Mediation models with moderate values of tan(/3) we find that 
the zero-lepton MET searches of supersymmetry at the LHC are only starting to be com- 
petitive with the Higgs bounds from LEP. From this perspective and for such models, the 
SUSY searches at the LHC are still very much in the beginning. 



1 Introduction 



The LHC has completed its first full year of collecting data and searching for new physics in 
7 TeV proton- proton collisions. The ATLAS and CMS collaborations conducted an impressive 
series of searches for supersymmetry (SUSY) at the LHC looking for squarks and gluinos in 
final states containing jets and missing transverse momentum [l}j4] . The latest of these analyses 
are based on ~ 1 fb~^ of data taken in the first half of 2011 |ln4r 

So far, these searches have shown no evidence for SUSY (nor for any other new physics). 
This fact by itself does not imply the demise of supersymmetry, it rather points out that new 
physics was not 'hiding around the corner' and more work is needed to uncover it. In this paper 
we will assess the impact of these searches for General Gauge Mediated (GGM) SUSY breaking 
models with long-lived next-to-lightest supersymmetric particles (NLSPs). Interpretations of 
the latest LHC data in terms of other phenomenological models are given in (6}{ll]. Here, 



we will derive the exclusion contours arising from the ATLAS zero-lepton search for SUSY at 
1.04 fb~^ [l], henceforth 'jets + MET', and compare with the LEP Higgs bounds fsj. We wih 
argue that for GGM models with low-to- moderate values of tan(/3), the Higgs bounds are still 
the strongest bounds in sizable areas of parameter space. As such for this classes of models 
supersymmetry searches at the LHC are only beginning. 



2 Main Results: Summary and Interpretation 

Figure [T] summarises our main findings for the LHC SUSY and LEP Higgs exclusion contours 
in models of General Gauge Mediationj^ in the context of the minimal supersymmetric standard 
model. We plot both sets of exclusions on the plane of physical squark (vertical axis) and 
gluino masses (horizontal axis). The left panel shows GGM with lower messenger masses of 
M = 10^ GeV, and the right panel shows a high-messenger-mass case with M = 10^^ GeV. 
The SUSY exclusion contours in GGM are derived from 1.04 fb~^ of data collected by the 
ATLAS collaboration in their searches for jets and missing transverse momentum (MET) in 
proton-proton collisions at 7 TeV with zero leptons in the final state [Ij. We obtain these SUSY 
exclusion contours following the same strategy as in (isl. As will be explained in more detail 



in the following section, we use Herwig++ 14,151 to simulate events in GGM and RIVET 16 



to implement the experimental cuts used by ATLAS in fl]. The resulting signal rates of GGM 
are then compared to data. 

The data strongly constrains events with missing energy carried for example by a neutralino 
NLSP. In sizable regions of GGM parameter space the NLSP is the stau (rather than the neu- 
tralino), which we take to be stable on collider scales. In our analysis we use a conservative 
approach to the data by assuming that events with long-lived stau NLSPs will not give MET 
signatures and thus will not be constrained. In these regions one should instead use the dedi- 



cated searches for long-lived charged particles. These searches have already started 17,18 and 
we plan to return to these analysis in forthcoming work. 

The red and blue solid lines in each panel of Fig. [T] denote these MET exclusions in GGM 
for low and for high values of tan(/3), with tan(/3) = 5 blue and tan(/3) = 45 given by the red 
curve. Following the contour from the upper left, it can be easily seen that the bounds from jets 
-|- MET searches are largely tan(/3)-insensitive until one reaches the turning point where the 



^We refer the reader to Sect. 3.3 and references therein for details on GGM. 
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Figure 1: Exclusion contours for General Gauge Mediation from the ATLAS 1.04 fb^^ jets 
+ MET SUSY search versus Higgs bounds from LEP in terms of physical squark and gluino 
masses. The left panel shows GGM with a messenger scale M = 10^ GeV, the right panel 
M = IQ^^ GeV. Jets + MET exclusions are shown as solid lines with tan(/3) = 5 in red and 
tan(/3) = 45 in blue. Higgs exclusion contours obtained from HiggsBounds are shown as dashed 
lines in red, black and blue, for tan(/3) = 5, 10 and 45 respectively. For tan(/3) = 45 there is a 
region where the stau becomes tachyonic, indicated by the dotted blue line. The grey area is 
theoretically inaccessible 12 



NLSP changes from the lightest neutralino to the stau. This point itself is tan(/3) dependent 
as can be seen from the figure. Beyond this point the curve follows the boundary between 
neutralino and stau NLSP regions. As can be seen from Fig. [l] the stau NLSP region grows 
at large tan(/3) because the stau mass is reduced by its Yukawa coupling, which increases with 
tan(/3). 

GGM models also allow for NLSPs decaying inside the detector. This can be either prompt 
or with a displaced vertex. This case has been studied, e.g. in An example is a bino-like 
neutralino decaying into a photon and gravitino. In this case it is more appropriate to rely 



on the dedicated MET searches with two photons in the final state 19 , 20 . These are more 



sensitive than the inclusive jets + MET searches, since the missing transverse momentum is 
halved when each neutralino can decay into a photon and a gravitino (with only the latter 
carrying the MET). In the present paper we concentrate on long-lived NLSPs. 
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Figure 2: Comparison of the ATLAS 1.04 fb~ jets + MET exclusion (solid line) with the LEP 
Higgs bounds (dashed hues) for the CMSSM. In the left panel we set = and show the 
Higgs bounds for tan(/3) = 5 in red, tan(/3) = 10 in black and tan(/3) = 45 in blue. On the right 
panel we fix tan(/3) = 10 and show Aq = —mo, 0, +2mo in green, purple and brown. 



The Higgs-based exclusion contours for GGM are obtained by using FeynHiggs 21-24 and 



HiggsBounds 25,26] (for details see the next section). This automatically enforces the most- 
constraining-to-date bound on the Higgs sector of each model. In our case this reduces to the 
LEP searches for SM-like Higgs bosons, so we will refer to these exclusions collectively as the 
LEP Higgs bounds. 

The LEP Higgs bounds for GGM are shown in Fig. [T] as dashed lines in red, black and 
blue, for tan(/3) = 5, 10 and 45 respectively. The position of the Higgs exclusion contours in 
Fig. [T] shows significant dependence on the value of tan(/3). This is as expected, since the tree 
level contribution to the mass of the lightest Higgs is m^'if'^ = mz cos(2/3), which decreases with 
decreasing tan(/3), making it harder to satisfy the LEP lower bound on the Higgs mass. In 
order to fulfill this bound, as is well known, one needs large loop corrections from the stop. At 



Quite clearly, to talk about Higgs exclusions contours one has to use non-trivial information about the SUSY 
spectra. This data is usually not contained in simplified models. Therefore it is desirable to use complete models, 
as we do here. 
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one- loop level (see e.g. [27]) one has, 
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(m| - m| ) ln(m| /m| ) } /m? 



,(2.1) 



where and are the cosine and sine of the stop mixing angle. 

To compensate for the tree level deficit at lower values of tan(/3) we need sufficiently large 
stop masses which makes the Higgs exclusion contours more powerful, i.e. they exclude larger 
portions of the parameter space. 

From Fig. [l] we see that in the GGM SUSY framework with not too high values of tan(/3), the 
current LHC searches for supersymmetry based on missing transverse momentum with jets and 
zero leptons in the final state have only just started to be competitive with the Higgs exclusions 
from LEP. In these scenarios we are not surprised that SUSY has not yet been discovered at 
LHC. 

We note however, that the jets + MET searches exclude even decoupled squarks for suf- 
ficiently light gluinos [rrig < 500 GeV). This is a region where the Higgs exclusion cannot 
compete. 

The area of SUSY parameter space least probed and constrained by jets+MET SUSY 
searches is the area with a charged NLSP, in GGM notably the stau. This has already been 
noted in 13 and can be seen very clearly from Fig. [TJ The Higgs exclusion holds regardless of 
the NLSP identity. Therefore, it is crucial that the data from dedicated searches for charged 



massive particles started in 17, 18 can be interpreted in terms of general classes of SUSY 
models. 



To what extent is the relative prominence of the Higgs exclusion in GGM different from 
other SUSY models, in particular from the much studied example of the CMSSM? An important 
distinctive feature of all gauge mediation models in this context is the absence (or more precisely 
the loop-level suppression) of the trilinear soft j4-terms at the high scale. This is not the case 
in gravity mediation models. In the CMSSM the A-coupling at the GUT scale is given by 
which is a free parameter. 

We now switch to the CMSSM and compare the jets + MET exclusion with the LEP Higgs 
exclusion at different tan(/3) and Aq. The insensitivity of the former to these variables in the 
CMSSM was shown in |28]. Higgs bounds however depend on both. We illustrate this in 
Figure [2| In the left panel we show tan(/3) = 5, 10 and 45 for constant Aq = 0, and we reach 
similar conclusions as for the gauge mediation case discussed above. In the right panel we 
vary Aq from —mo to +2mQ (where rriQ is the usual CMSSM scalar mass parameter), keeping 
tan(/3) = 10 fixed. The variation of Aq affects the value of the Higgs mass and the relevance 
of the Higgs exclusion contour. The reason is that negative (positive) Aq increases (decreases) 



stop mixing and thus the physical Higgs mass (see Eq. (2.1)), making the Higgs searches less 
(more) constraining. 

In this sense, for the CMSSM and related models, the LEP Higgs bound provides a stronger 
constraint than jets + MET searches only in a part of the parameter space, characterised 
by vanishing or (moderately) positive values of ^o- A common strategy for dealing with the 
additional parameter Aq is, instead of keeping it fixed, to optimise with respect to it. In other 
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Figure 3: Jets + MET SUSY exclusion for the CMSSM, and for GGM for different values 
of tan(/3) and the messenger scale M. The diagonals delimiting the excluded areas in each 
case arise from one of the following two effects. Either the NLSP (in gauge mediation models) 
changes from a neutralino to a stau, for which our implementation of the ATLAS jets + MET 
search has no sensitivity, or we simply reach a region in the gluino-squark mass plane which is 
theoretically inaccessible (12 



6 



words one looks for such values of Aq where the SUSY is least constrained by the Higgs (and, 
if desired, other experimental inputs [7||29]). For GGM models, Aq is not a free parameter and 
there is no simple way of reducing the dominance of the LEP Higgs exclusion contours (apart 
from dialing up tan(/3)). 

We would like to note that even in models where Aq is a free parameter, the direct searches 
provide additional sensitivity beyond the Higgs bounds only for a certain range of Aq. As the 
overall scale probed by the collider increases this range grows as well. Therefore the size of the 
newly probed parameter space increases rapidly. Again SUSY searches are just taking off. 

The observation that for not too high values of tan(/3) < 10, the available parameter space 
of general models of gauge mediation is not yet significantly constrained by the current MET 
SUSY searches at the LHC is one of the main points of this paper. The ongoing searches of 
supersymmetry at the LHC are entering the crucial discovery phase, but they are far from 
announcing the demise of supersymmetry. 

The details of our analysis leading to the implementation and interpretation of the LHC 
SUSY searches for GGM models as well as the validation of this implementation for the known 
CMSSM case will be presented in the following section. There we also explain how the Higgs 
bounds were derived. Before moving on, however, we note that the jets + MET exclusion shows 
little model dependence provided the neutralino is stable on detector scales. In Figure [3] we 
show this exclusion in terms of physical squark and gluino masses both for the CMSSM and 
for GGM with low and high values of tan(/3) and of the messenger scale M. As one can clearly 
see the variations in the relative positions of the contours are not dramatic]^ The position of 
the turning point where stau becomes lighter than the neutralino and the resulting slope do 
however depend on the model. 



3 Details of the Analysis 



3.1 Implementation and validation of the ATLAS jets + MET SUSY search 

Our approach in implementing and validating the latest 1.04 fb^^ ATLAS MET search con- 
straints for a general BSM model follows the general strategy laid out in Ref. 13 and uses a 
combination of Herwig++, RIVET and Prospino. 

A Monte Carlo event generator, in our case Herwig++ 2.5.1 [l4|[T5], is used to compute 
signal rates from any given BSM model, in our case GGM (or CMSSM for validation) with 
mass spectra calculated by SoftSUSY 3.1.6 ^31j. Our implementatiorj^ of the experimental event 
selection used by ATLAS in [l], defining five different signal regions, was done in the RIVET 
1.5.2 [l6] analysis framework. Since the matrix elements implemented in general purpose Monte 
Carlo event generators are only accurate to leading order in perturbative QCD, we supplement 
the acceptance calculated with Herwig++ with an NLO cross section calculated with Prospino 
2.1 32-35 . More precisely, we use Herwig++ to simulate three sets of supersymmetric particle 



production processes for each point in SUSY parameter space: a) squark and gluino production, 
b) the production of an electroweak gaugino in association with a squark or gluino and c) 
the production of slepton and electroweak gaugino pairs. The fraction of events passing the 



•^This agrees with the observation 
''it is pubhcly available as part of 



30 that the slepton mass has little effect on the exclusion. 
RIVET package: analysis ATLAS.2011.S9212183. 
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Figure 4: Validation of our signal simulation in the CMSSM against the ATLAS limits ob- 
tained in [l] in the five signal regions defined in their Table 2: > 2-jets, > 3-jets, > 4-jets with 
rUes > 500 GeV, > 4-jets with rries > 1000 GeV and "High Mass" counting from left to right 
and top to bottom. We use an efficiency factor as described in the text. The red curves denote 
the exclusion limits obtained from our simulation; its scale variation (changing the renormal- 
ization/factorization scale by 2^^) is indicated by the black curves. The ATLAS exclusion is 
shown in yellow. In the last panel - the High Mass region - we also show the ztlcr contours for 
the ATLAS expected exclusion (thin dashed white) . The wide splitting of these two contours 
shows the difficulty in predicting this contour's position precisely. The colour indicates the cross 
section in fb in each of the signal regions. The red line in the colour bar shows the quoted limit. 
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experimental cuts in each channel was then used together with the cross section calculated 
using Prospino to obtain the number of signal events passing the cuts for each of the five signal 
regions. The maximum allowed number of signal events in each signal region is provided by 
ATLAS [l]. Comparing our numbers of events calculated in a given SUSY model point with 
this number decides whether the model point is allowed or excluded. 

To validate our implementation we simulate the same CMSSM plane in the five signal regions 
presented by ATLAS in Ref. [l] and compare our exclusion with theirs. As can be seen in Fig. |4] 
we obtain excellent agreement between our exclusion contours and theirs. Having implemented 
the same kinematic cuts and selections as ATLAS this consistency means that the detector 
efficiencies are relatively close to one in the relevant kinematic regions. This makes us confident 
that our analysis, which does not use a detector simulation, can be applied to the GGM case. 
As discussed in Ref. [l| an efficiency factor is needed to account for loss of jet energy. We have 
included a factor of 0.85 for the four of five signal regions with softer cuts on the jet px (i.e. 
excluding the High Mass regiorr| to account for this and give a minor improvement. 



3.2 Higgs Bounds 

To check exclusion by Higgs-based searches, we used the superpartner spectrum computed by 



SoftSUSY with FeynHiggs 2.8.5 121-24 to accurately calculate further details of the Higgs sector, 



in particular decay cross-sections. This information was passed to HiggsBounds 3.5. Obeta [25 



26 to compare with current experimental limits from LEP, the Tevatron and the LHC (ver- 
sion 3.5.0 includes data from the EPS-HEP-2011, LeptonPhoton-2011, and SUSY-2011 con- 
ferences). HiggsBounds returns 0"/<Jiimit for the search channel with the highest statistical 
sensitivity: a ratio greater than 1 indicates 95% confidence-level exclusion. This approach 
has the rigour of Monte Carlo based exclusion - in comparing cross-sections to their limits 
rather than using model-dependent mass limits - whilst only taking a time comparable to 
the generation of the mass spectrum itself. In our case we see only the LEP lower bound, 
m/j > 114.4 GeV, for an SM-like Higgs. However for more exotic SUSY scenarios, and with 
more data collected by the LHC and included in HiggsBounds, this may not be the case. (A 
user- friendly python script for linking SoftSUSY with FeynHiggs with HiggsBounds is available 
at [http: / / www.ippp.dur.ac.uk/^SUSY, ) 



3.3 Implementation of the ATLAS jets + MET SUSY search in GGM 



In General Gauge Mediation [36] the soft supersymmetry breaking gaugino masses at the mes- 
senger scale M are given by 

M^^(M)= k^^^Kc^i (3.1) 

where ki = (5/3, 1, 1), kiai (no sum) are equal at the GUT scale and are the gauge coupling 
constants. Similarly, the scalar mass squareds are 

,.2 1 



m){M) = 2Y,C.h^(^ (3-2) 



where the Cj are the quadratic Casimir operators of the gauge groups. 

^For the High Mass region our agreement with the ATLAS exclusion contour is better without such an 
efficiency factor. As can be seen from Fig. |4] the position of this contour has large uncertainties. 
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Figure 5: Exclusion contours for General Gauge Mediation from the ATLAS 1.04 fb^ jets + 
MET SUSY search and from the LEP Higgs bound, with a messenger scale M = lO'^ GeV (top 
panels) and M = 10^^ GeV (bottoms panels). The left panels have tan(/3) = 5 and the right 
panels tan(/3) = 45. The red line shows implementation of the ATLAS data for GGM; the 
black lines denote the scale variation. The yellow dashed line is the Higgs exclusion. White 
regions are stau NLSP. The blue dotted line in the bottom right panel delineates a tachyonic 
stau region. The colour throughout each plot denotes cr/aw^ait for the most constraining region. 
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The value of the messenger scale M, the value of tan(/3) together with Aq and A5 appearing 



in Eqs. (3.1)-(3.2) at M characterise a point in the GGM parameter space [37| . 

For clarity we show numerical results for the case where, as in [38','39] , we do not further split 
Ag or As into separate independent factors for each gauge group (this also ensures automatic 
gauge unification). At the end of this section we explain, however, that our findings hold 
also for large regions of the full GGM parameter space with separate A for each gauge group. 



Furthermore, and in contrast to [38,39 , we do not confine ourselves here to pure GGM models 
where tan(/3) is determined from the vanishing input value of the B soft term. Here tan(/3) is 
treated as a free parameter. 

Before we proceed let us briefly comment on the nature of the NLSP in gauge mediatioiij^ 
As is well known the LSP in gauge mediation is the gravitino. The collider phenomenology 
then depends crucially on the nature of the NLSP and its lifetime. In general gauge mediation 
there is a wide range of options for what can be the NLSP. However the jets + MET search 
as implemented in this paper is most sensitive to long-lived neutralino NLSPs. A long lived 
neutralino naturally leads to significant amounts of missing energy. On the contrary a non- 
neutral NLSP like, e.g., a stau would leave a track and therefore wouldn't be counted as missing 
energy. Similarly if the neutralino decays inside the detector into, say, a photon and a gravitino, 
then less missing energy is generated (and the event may also be vetoed due to the presence 
of photons). For these reasons the zero-lepton jets + MET searches are primarily aimed at 
and are most sensitive to long-lived neutralinos. This type of NLSP arises in large parts of the 
GGM parameter space. Specifically, the NLSP decay length is approximately given by, 

1 /lOOGeVy/^A^V., 
kQ V "T-NLSP / VlOOTeVy 

where kc ■= 1/Ggrav quantifies the coupling of messengers to the SUSY breaking sector which 
can take values of order 1 but is often much smaller. For Ac ~ 10^ GeV, already moderate 
messenger scales of the order of lO'' GeV, lead to Ldecay ^ 10 m and a decay outside of the 
detector. When /cc ^ 1 even smaller values of the M may suffice. This addresses the longevity 
of the NLSP. Turning now to the NLSP species we find in the simple setup with a single Aq 
and a single A^ that roughly half of the parameter space is neutralino (the rest is mostly stau). 
The NLSP content of the current model is very similar (especially at high values of tan/3) to 
that of pure GGM which is shown in Figure 5 of Ref . [39] . 

Figure [5] shows our results for the jets + MET exclusion of GGM in terms of the original 
Ac and A5 model parameters. The upper panels show a messenger scale M = 10^ GeV, the 
lower panels M = 10^^ GeV. The left panels have tan(/3) = 5 and the right panels tan(/3) = 45. 

The relative importance of Higgs (yellow dashed line) vs direct SUSY exclusion contours in 
GGM for various values of tan(/3) (and M) can be directly read off these Aq-As planes. The 
projection of Fig. [5] to the physical mg-msq plane is Fig. [ij 

Let us conclude with a brief discussion of what happens if we allow for six different As 
to span the entire GGM parameter space. We need to address what happens in this case for 
both the direct SUSY search and Higgs exclusion contours. Jet + MET SUSY searches depend 
crucially on three parameters, the mass of the gluinos, the first generation squark mass and 



®A detailed study for the simpler case of pure GGM can be found in Sect. 3 of Ref. 
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the neutralino mass. The first two are taken into account by presenting our results in Fig. [T] 
on the gluino-squark mass plane. The gluino mass chiefly depends on the parameter Ag,3- For 
direct searches the most relevant squark mass is the lightest one. This squark mass squared is 
determined by a certain linear combination of the and the ■ through the RG evolution. 
As long as the neutralino remains the NLSP and its mass is not fine-tuned to be very close 
to the squark and gluino masses, the exclusion limits from direct searches in the gluino-squark 
mass plane are not affected by a splitting of As. On the other hand if the splitting of As changes 
NLSP species jets + MET searches lose sensitivity, as explained earlier. 

Now the Higgs exclusion contour chiefly depends on the stop masses. Since gauge mediation 
does not distinguish between the generations this mass is directly linked to the mass of the first 
generation squarks. Unless extreme splittings are introduced this is essentially the same squark 
mass parameter that governs the direct searches. In addition there is a higher order dependence 
of the Higgs mass on the gluino mass, which is our second parameter. In summary the same 
two parameters, squark and gluino mass, are the most relevant ones for both, the direct SUSY 
searches and the Higss constraints. Showing the results in the plane spanned by these two 
parameters gives qualitatively the same picture independent of the particular combination of 
the original high energy parameters, A. 

Note added 

To bring the Higgs mass in our GGM model up into the region rrihQ ~ 125 GeV, hinted at by 
the recent LHC results ^40||41| , requires fairly high superpartner masses. If we want to keep 
gluino masses as low as possible (~ 500 GeV, cf. Fig. [T]), we find that the squark masses must 
be in the region of > 10 TeV. 
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